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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  In  this  report  can  be  converted  to 
Metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

Inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  Inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10“3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C) 

temperature  readings 

from  Fahrenheit  (F)  readings. 

use  formula:  C  -  (5/9) 

(F  -32). 

To  obtain  Kelvin  (K) 

readings,  use  formula: 

K  -  (5/9)  (F  -32)  +  273.15. 

SYMBOLS  AND  DEFINITIONS 

E  total  variance  in  wind  sea,  often  called  energy 
E(f)  variance  density,  often  called  energy  density 
E^  depth-limited  value  of  total  variance 

E  (f)  upper  bound  on  energy  density  in  a  frequency,  f 

n 

F  variance  density  spectrum  in  wave  number  space 

f  frequency 

f  low-frequency  cutoff 

c 

f  peak  frequency  of  the  spectrum 

H  depth-controlled  wave  height  (spectral) 

depth-limited  wave  height  (monochromatic) 

H 

t  depth-limited  wave  height  (irregular  sea) 

%o  zero-moment  wave  height,  also  called  significant  wave  height 
H  largest  individual  wave 

Hi/ 3  significant  wave  height 

h  depth 

k  wave  number 

R  transcendental  function  of  dimensionless  frequency 

U  windspeed 

a  Phillips'  equilibrium  coefficient 

ir  3.1415 

4  dimensionless  function  describing  deviation  from  deepwater  equilibrium 

range 

u.  dimensionless  combination  of  g,  f,  and  h 
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DEPTH-LIMITED  SIGNIFICANT  WAVE  HEIGHT:  A  SPECTRAL  APPROACH 

by 

C.  Linaood  Vincent 
I .  INTRODUCTION 

Research  into  the  shape  of  wind  wave  spectra  in  finite-depth  water  has 
suggested  an  expression  for  the  upper  limit  on  the  energy  density  as  a  func¬ 
tion  of  depth  and  frequency  (Kitaigorodskii,  Krasitskii,  and  Zaslavskii,  1975). 
In  this  report  this  expression  is  integrated  over  the  part  of  the  spectrum 
expected  to  contain  energy  to  estimate  a  limit  on  the  energy,  E,  in  the  wind 
wave  spectrum  and  to  define  a  depth-limited  significant  wave  height,  H£t  „ 

-  4.0(E)1/2  (1) 

More  precisely,  the  quantity  estimated  is  the  variance  of  the  sea  surface  to 
which  E  is  directly  related.  Following  convention,  E  and  E(f)  denote 
energy  and  energy  density  spectrum  although  the  true  units  of  computation  are 
length  squared  and  length  squared  per  hertz.  The  term  zero-moment  wave  height, 
%o  will  be  used  to  denote  4.0(E)1/2.  H1/3  is  the  average  height  of  the  one- 

third  highest  waves.  denotes  values  of  HMq  that  are  depth  limited.  In 

deep  water,  H^0  is  approximately  H1/3,  but  this  is  not  necessarily  true  in 
shallow  depths.  refers  to  the  depth-limited  monochromatic  wave.  The  vari¬ 

ation  of  H/  with  depth,  h,  is  Investigated  and  compared  with  the  mono- 
chromatically  derived  depth-limited  wave  height,  H^.  Because  HMq  and  Hi/ 3 
are  about  equal  in  deep  water,  they  are  both  frequently  called  significant 
wave  height. 

This  report  briefly  reviews  the  theoretical  development  of  the  limiting 
form  for  spectral  densities  as  a  function  of  water  depth  and  presents  field 
evidence  supporting  this  form.  The  simple  derivation  of  the  depth-limited 
energy  and  significant  wave  height  is  then  given,  followed  by  field  and  labora¬ 
tory  data  evaluating  the  prediction  equation.  Unless  otherwise  noted,  the 
developments  of  this  report  are  restricted  to  wave  conditions  described  by  a 
wave  spectrum  of  some  width  such  as  an  active  wind  sea  or  a  decaying  sea. 

II .  THEORETICAL  BACKGROUND 

Phillips  (1958)  suggested  that  there  should  be  a  region  of  the  spectrum 
of  wind-generated  gravity  waves  in  which  the  energy  is  limited  by  wave  steep¬ 
ness.  Phillips  derived  an  expression  for  the  limiting  density  in  deep  water: 

Ea(f)  -  ag2f"5(2irr4  (2) 

where  a  was  considered  to  be  a  universal  constant.  Field  studies  reviewed 
by  Plant  (1980)  demonstrated  that  equation  (2)  adequately  describes  the  part 
of  the  wind  sea  spectrum  above  the  peak  frequency  of  the  spectrum.  However, 
Hasaelmann,  et  al.  (1973)  indicated  that  the  equilibrium  coefficient  a  is 
not  constant  but  varies  systematically  with  wave  growth  leading  the  authors 
to  speculate  that  resonant  interactions  in  the  spectrum  force  the  spectrum  to 
evolve  to  the  form  of  equation  (2).  Toba  (1973)  suggested  that  the  equilibrium 
range  form  might  be  proportional  to  U*f-t*  in  order  to  remove  the  variation 
of  a. 
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Kitaigorodskii ,  Krasitskii,  and  Zaslavskii  (1975),  using  Phillips'  (1958) 
expression  for  the  steepness  limited  form  of  a  wave  spectrum,  F,  in  terms  of 
the  wave  number  modulus, 


F(k) 


r”3 


(3) 


solved  the  transformation  of  equation  (3)  to  a  frequency  spectrum  in  finite- 
depth  water.  The  finite-depth  form,  E  (f,h),  was  shown  to  be  equal  to  the 
deepwater  form  (eq.  2)  times  a  dimensionless  function,  4(si^), 


Em(f*h)  "  ff^  ‘K* 


(4) 


Kitaigorodskii,  Krasitskii,  and  Zaslavskii  suggested  a  value  of  0.0081  for  a. 

The  function  4  requires  an  iterative  procedure  for  solution  and  is 
defined  as 


*(u>h)  -  R-2  (o^)  I  1  +  — 


[' 


with 
where  oi 


Sinh  (2u^)R(«^)J 

«h  -  ^(h/g)1/2 

2irf  and  R(w^)  is  obtained  from  the  solution  of 


Rirn^)  tanh^u2  R^)^ 


(5) 

(6) 

(7) 


The  dimensionless  parameter  u.  related  the  frequency  and  depth  to  the  devia¬ 
tion  from  the  deepwater  form.  "Whep  is  greater  than  2.5,  *  is  approxi¬ 
mately  1,  when  tit^  is  zero  *  is  zero;  When  is  less  than  1 

♦(w^)  as  Wjj/2  (8) 

For  less  than  l,  a  combination  of  equations  (8)  and  (4)  leads  to  the 

expression 

E  (f ,h)  -  aghf"3/(2(2ir)2)  (9) 

m 

Thus  in  the  shallow-water  limit,  the  bound  on  energy  density  in  the  wave 
spectrum  is  proportional  to  f”3  compared  to  f~5  in  deep  water,  and  depth 
is  included  linearly. 

Resio  and  Tracy  (U.S.  Army  Engineer  Waterways  Experiment  Station,  personal 
communication,  1981)  have  analyzed  the  resonant  interactions  and  derived 
equivalent  expressions  to  equations  (3)  and  (4)  on  the  basis  of  similarity 
theory.  The  conclusion  of  their  theoretical  study  is  that  the  role  of  the 
wave-wave  interactions  in  both  deep  and  shallow  water  is  to  force  the  spectrum 
to  evolve  to  the  form  of  equation  (4).  Their  theory  may  be  distinguished  from 


that  of  Kitaigorodskii ,  Krasitskii,  and  Zaslavskii  (1975)  in  that  their 
coefficient  a  is  expected  to  vary  with  vave  conditions  and  not  remain  a 
universal  constant. 

III.  FIELD  EVIDENCE  FOR  THE  FINITE -DEPTH  SPECTRAL  FORM 

Prior  to  Kitaigorodskii,  Krasitskii,  and  Zaslavskii  (1975),  Kakimuma 
(1967)  and  Druat,  Massel,  and  Zeldler  (1969)  had  noted  that  the  shape  of  the 
spectrum  in  shallow  water  deviated  'from  Phillips'  (1958)  form.  Kitaigorodskii, 
Krasitskii,  and  Zaslavskii  cited  evidence  supporting  the  f "3  form,  as  did 
Thornton  (1977)  and  Gadzhiyev  and  Kratsitsky  (1978).  Ou  (1980)  provided 
laboratory  evidence  for  equation  (4) .  A  review  of  spectra  collected  at  the 
Coastal  Engineering  Research  Center's  (CERC)  Field  Research  Facility  (FRF)  at 
Duck,  North  Carolina,  and  at  other  gages  in  shallow  water  supports  a  near 
f-3  spectral  slope  in  depths  less  than  10  meters  for  large  wave  energies. 

These  findings  indicate  a  further  evaluation  is  needed  of  how  well  the 
equation  fits  observed  spectra.  During  the  Atlantic  Remote  Sensing  Land  and 
Ocean  Experiment  (ARSLOE)  conducted  in  October  and  November  1980  at  the  FRF, 
North  Carolina,  wave  spectra  were  collected  in  36  meters  of  water  about  36 
kilometers  offshore  of  the  CERC  facility  (Fig.  1),  using  the  National  Ocean 
Survey's  directional  buoy,  XERB,  with  accelerometer  buoys  in  depths  of  25,  18, 
and  17  meters  of  water  located  at  distances  of  12,  6,  and  3  kilometers  offshore 
along  a  line  from  the  facility  to  the  XERB.  In  addition,  data  from  Baylor 
gages  at  seven  locations  in  1.5-  to  9-meter  depths  along  the  FRF  pier  were 
collected.  On  25  October  1980  a  large,  low-pressure  system  generated  waves 
with  significant  heights  up  to  5.0  meters.  Data  were  collected  continuously 
at  the  XERB.  during  the  per ’od  of  high  waves  and  spectra  at  all  gages  were 
computed  every  20  minut.  >. 

As  a  test  the  obsr  -ved  spectra,  E(f),  were  normalized  to  the  following 

forms 

a  (f)  =  f5E(f)  (2irr4/g2  (10) 

P 

(f)  -  f®E(f)  (2ir)4/g2  *(0^)  (11) 

a3  (f)  -  f3E(f)/gh  (12) 

Equation  (10)  is  an  estimate  of  the  equilibrium  coefficient  as  a  function  of 
frequency  if  the  spectra  follow  the  deepwater  form.  Likewise,  equation  (11) 
is  an  estimate  of  the  coefficient  if  the  spectra  follow  the  proposed  finite- 
depth  form,  and  equation  (12)  is  an  estimate  of  the  coefficient  if  the  proposed 
shallow  water  (u^  less  than  1)  holds  over  most  of  the  spectrun.  If  any  of 
these  forms  fit  aspectrum  then  the  corresponding  function  a(f)  should  be 
constant  with  frequency.  Therefore  in  a  regression  of  f  against  a(f),  f 
should  explain  no  variance;  consequently,  the  degree  of  fit  to  the  spectrum 
by  each  of  the  three  forms  can  be  estimated  by  how  poorly  f  explains  variance 
in  the  regression  and  how  flat  the  slope  with  f  is.  The  regressions  were 
performed  over  the  region  from  the  spectral  peak  to  twice  the  spectral  peak 
and  the  results  are  tabulated  in  Tables  1  and  2. 


Table  1.  Normalized  form  regression  (analysis  average 
of  percent  variance1  in  regression  of 
normalized  form  against  frequency,  f, 
explained  by  f). 


Gage 

Mean 

depth 

(m) 

Deepwater 

form2 

Finite-depth 

form3 

Shallow-water 
limit  form4 

XERB 

36 

43 

17 

48 

710 

25 

53 

28 

8 

630 

18 

40 

13 

10 

620 

17 

55 

24 

18 

625 

68 

35 

32 

655 

70 

60 

59 

615 

72 

60 

59 

1Since  the  proposed  form  is  supposed  to  remove  variation 
with  f,  a  high  explained  variance  with  f  indicates  that 
the  form  does  not  fit  the  spectra  well. 


2f5  E( f)  (2ir)‘4/g2 

3f5  E(f)  (2ir)“4/g2  4  (<*) 

4f3  E(f)/gh 


Table  2.  Average  slope,  X10"3, 
against  f 1 . 


Gage 

ESI 

Oh<f) 

03(f) 

XERB 

31 

17 

10 

710 

54 

35 

4 

630 

47 

18 

-3 

620 

49 

25 

3 

625 

82 

81 

35 

655 

60 

101 

53 

615 

18 

67 

33 

1 Slope  a  is  unit  change  of 
a  per  hertz.  Region  of  the 
spectnss  analyzed  in  regression 
analysis  is  about  0.1  hertz. 


Using  the  criteria  established  above,  data  summarized  in  Table  1  indicate 
that  in  all  cases  either  the  finite-depth  form  or  the  f-3  limit  appears  to 
fit  the  spectra  better  than  the  deepwater  form.  This  is  because  f  consist¬ 
ently  explains  less  variance  in  these  regressions  than  in  the  regressions 
against  the  deepwater  form.  In  a  regression  analysis  under  an  assumption  of 
normally  distributed  variates,  the  hypothesis  of  zero  correlation  is  rejected 
for  the  number  of  frequency  components  from  fp  to  2fp  if  the  regression 
coefficient  is  greater  than  0.632  at  a  5  percent  level  of  significance.  This 
translates  to  a  value  of  40  percent  for  the  values  in  Table  2.  Table  1  indi¬ 
cates  that  the  average  R2  for  the  regressions  in  the  deepwater  form  are 
always  greater  than  40  percent,  suggesting  that  there  is  correlation  with  f . 
The  average  finite-depth  form  value  is  less  than  40  percent  for  all  but  two 
(655  and  615)  of  the  gages,  suggesting  a  tendency  for  no  correlation  with  f. 
The  shallow-water  limit  results  suggest  zero  correlation  except  for  gages  XERB, 
655,  and  615.  Table  2  indicates  that  the  slopes  are,  in  general,  lower  as 
well.  Plots  of  f3E(f)  and  fsE(f)  show  that  the  spectra  appear  to  more 
closely  follow  a  f“3  slope  (Fig.  2). 

The  results  of  the  regression  analysis  for  the  gages  at  depths  greater 
than  9  meters  appear  to  be  more  closely  fit  by  a  f"3  form  than  the  results 
at  9  meters  and  at  shallower  gages.  The  observed  spectra  at  the  shallower 
gages  tend  to  be  less  than  the  proposed  upper  limit.  It  is  thought  that 
refraction,  bottom  friction,  and  massive  breaking  must  dominate  the  spectra 
in  and  around  the  peak,  suppressing  the  values  below  the  proposed  limiting 
value.  This  would  indicate  that  in  very  shallow  water,  the  proposed  form  may 
be  conservative.  Plots  of  storm  spectra  at  different  gage  sites  are  compared 
to  the  limiting  form  in  Figure  3. 

The  variation  of  the  equilibrium  coefficient  ot  computed  over  the  range 
fp  to  2fp  varies  based  on  gage  and  time  (as  represented  by  sea  and  swell 
conditions),  with  a  for  the  sea  conditions  being  larger.  Additionally, 
there  appeared  to  be  a  tendency  for  a  to  increase  slightly  from  deep  to 
shallow  water.  On  occasion  a  calculated  at  the  peak  of  the  spectrum  exceeded 
the  value  of  0.0081.  However,  when' the  a  value  at  the  peak  was  compared  to 
the  a  value  averaged  over  the  frequencies  from  fp  to  2f ,  it  was  evident 
that  the  average  value  was  much  less  than  the  value  at  the  peak. 

The  field  evidence  from  a  variety  of  sources  supports  the  conclusion  that 
the  maximum  energy  densities  above  the  peak  frequency  of  the  spectrum  can  be 
approximated  by  equation  (4),  which  in  the  shallow-water  limit  approaches 
equation  (9).  Evidence  from  Ou  (1980)  and  the  data  in  this  report  suggest 
that  the  coefficient  a  may  not  be  a  universal  constant.  There  is  also 
evidence  that  once  very  shallow  depths  are  reached,  other  mechanisms  can 
dominate  spectral  shape  in  the  vicinity  of  the  peak;  the  deviation,  however, 
is  such  that  equation  (4)  appears  to  be  an  overestimate. 

IV.  FORMULATION  OF  DEPTH-LIMITED  SIGNIFICANT  WAVE  HEIGHT,  H£ 

Since  equation  (4)  provides  an  estimate  of  the  upper  limit  on  energy 
density  in  water  depth  h  as  a  function  of  frequency  and  wave  generation 
condition  as  expressed  by  the  coefficient  a,  it  is  possible  to  estimate  the 
upper  bound  on  the  depth-limited  wave  energy,  E^,  if  a  low-frequency  cutoff 
value,  fc,  is  known.  Ej^  can  simply  be  estimated  by 


-3  Qep®n^®nC® 


E 

03 

W 

2 

8 

a 

O 

o  j 

(M 

• 

’e» 

1 

w  *1 

O  P* 

#y  A  1 
O..M  * 

ro 

u3 

*0  I 

M-4 

®  * 

4 

+ 

\r> 

ta 

$2 

r> 

'»  % 

1  -  s 

Uj 

II 


00-2* 


M*',S  O  H  » 

«  *0  £  5  « 

_  '*’*  O  fl>  "O 

|SS|«* 

s  S  »*  a 

.  .rt  O  t> 

tt  H  «£  4 

* 1  q)  £  *>*g 

o  s*£  u  & 

m  O  *  *  p  »H 

L  «  «  2  * 

n>  g  2  «  «  rl 

§  *  %*  rl 

-  I**4- * ; 

,  «.  o  **  ’'* 

Sg  t*  $ 

•S  o  fl>o  ®  “■* 
(j-°$  g  •  . 

•3  «r*  ^  8  ^ 

t4  <  £  o 
•H  .  2  d  H 
*J  r-»  M  &  g*i 

5  *  S  0*5  S 

«  gS  •  s  di 

2  ",  o  •  %.$  p* 

R  r  -  .j  u  •*. 


4) 

ob  » 
o*-* 

5  «s 

£ 

13 

Is 

Is 


^  d  JJ  d  IT*'1  • 

as 

g  44  S.  L  C  ^  ?  Hi  -H 

;£u  r  R  u  ^  h 


Q,  |4  c  ,w  «  .* 

&sis.r, 

>  "  a 41  jb  °  *• 

2  S'  n  ”5  « 

'S**  I  lh 

°  4J  \  „  £1  «* 

*  *  44  S  ® 

O  ™  CBJ3  , 

O  ^  r-»  c 

•S !  sT  -s2' 

|2»  “  11 
o  g  (d  -o  ** 


v  *2 

|8  #1 


_  -r— - TrS-  t  41  31  90  ’ 


Limit  on  Enorgy  Donoity  PropotoO  by 
Kitoigorodtkif,  Kraoit«kiit<m4  Zotiovokli  (19751 


(13) 


Eh  - 


E»(f,h)df 


Ctg2  f‘5 


*(uh)/(2*)4  df 


The  depth-limited  significant  wave  height  (spectral)  is  then 


*  4.0  (E^1/2  (14) 

In  shallow  water,  H#  is  expected  to  be  different  from  H1/3,  but  how  differ¬ 
ent  is  uncertain.  Although  Hj /%  has  a  long  tradition  of  use  in  coastal  engi¬ 
neering,  the  wave  height  H£  defined  in  equation  (14)  appears  to  be  a  more 
consistent  parameter  because  it  is  directly  related  to  the  energy  of  the  wave 
field. 

Figure  4  provides  curves  of  H£  as  a  function  of  cutoff  frequency,  fc, 
and  depth,  h,  for  a  *  0.0081.  If  a  is  different  an  estimate  of  H£  for 
that  a  can  be  made  by 

H t  -  H^a/O.OOSl)1/2  (15) 

where  H£  is  estimated  with  a  of  0.0081. 

Clearly  the  cutoff  frequency  and  the  value  of  a  are  crucial  for  obtaining 
estimates  of  H£.  An  examination  of  storm  spectra  indicates  that  the  spectral 
peak  is  quite  sharp.  Consequently,  a  reasonable  choice  for  fc  would  be  about 
90  percent  of  fp.  If  there  is  evidence  of  more  energy  on  the  forward  face  of 
the  spectrum,  fc  could  be  estimated  by  using  a  lower  percentage.  The  param¬ 
eter  a  can  be  obtained  by  fitting  equation  (4)  to  observed  data  if  available. 
For  field  engineers,  most  often  this  may  not  be  possible  in  which  case  a  can 
be  estimated  by  knowledge  of  the  peak  frequency,  fp,  and  windspeed,  U, 
through  the  relationships  developed  by  Hasselmaun,  et  al.  (1973).  The  values 
of  fp  and  U  can  be  obtained  from  hindcasts  or  measurements.  Figure  5 
provides  values  of  (a/0.0081) */2  as  a  function  of  fp  and  U. 

When  the  primary  frequency  components  containing  the  major  part  of  the 
energy  are  in  shallow  water,  as  determined  by  the  condition  uh  <  1,  then  Em 
is  given  by  equation  (9).  This  can  be  integrated  analytically  to  give  an 
estimate  of  H  for  a  ■  0.0081 

H*  -  i  (agh)1/2  fc“l  (16) 

Equation  (16)  has  the  remarkable  consequence  of  suggesting  that  H£  defined 
as  4.0(E)1/2  varies  with  the  square  root  of  depth  when  the  primary  spectral 
components  are  depth  limited.  The  monochromatic  depth-limited  wave  height, 

Hd,  varies  linearly  with  h. 

V.  FIELD  AND  LABORATORY  EVIDENCE  FOR  DEPTH-LIMITED  SIGNIFICANT  HEIGHT,  H £ 


In  order  to  test  the  applicability  of  equations  (15)  and  (16)  in  predicting 
in  shallow  water,  laboratory  data  taken  by  Seelig  and  Broderick  (1981)  in 
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Figure  4.  Depth-limited  significant  wave  height,  H£,  as  a  function 
of  water  depth  and  cutoff  frequency.  Curves  are  calculated 
for  a  ■  0.0081.  is  plotted  for  lower  limit  of  0.8  h. 


Peak  Frequency  ,  fp  (Hz) 

Plot  of  R  «  (a/0.0081) 1^2  as  function  of  peak  frequency 
of  spectrum  and  windspeed  U.  Data  based  on  JONSHAP  wind 
sea  relationships  (Hasselaann,  et  al. ,  1973).  Coefficient 
R  is  used  to  adjust  curves  in  Figure  4  to  account  for 
variation  in  a . 


a  flume  44  meter*  long  and  0.45  meter  vide  with  a  maxima  water  depth  of  0.6 
meter  and  a  bottom  slope  of  1:30  at  one  end  were  examined.  Seelig  and  Broderick 
ran  a  variety  of  spectral  shapes  and  energies.  Figure  6  Is  a  plot  of  H» 
calculated,  as  in  equation  (1),  from  a  Fourier  analysis  of  their  wave  data 
against  h1^2.  Typically,  the  wave  appear  to  shoal  with  decreasing  depth, 
thereby  Increasing  in  height  until  a  point  is  reached  at  which  the  wave  height 
decreases  linearly  with  the  square  root  of  depth.  Figure  7  is  an  estimate  of 
H^,  based  on  equation  (15),  for  two  fonts  of  fc.  A  plot  of  the  maximum 
individual  wave,  H^g,  is  plotted  as  is  the  monochromatic  breaking  limit  which 
Hoax  appears  to  follow.  H£  is  much  less  than  the  monochromatic  breaking 
limit  in  this  case.  Figure  8  provides  plots  of  versus  h1/2  for  wave  data 
at  FRF  on  25  October  1980.  The  value  of  h  is  estimated  by  an  average  of 
profiles  before  and  after  the  storm  and  Includes  the  tide  and  the  wave  setup. 

The  curves  are  approximately  linear  with  h1'2. 

VI .  DISCUSSION 

An  examination  of  the  characteristics  of  spectral  shape  in  shallow  water 
has  led  to  a  method  of  estimating  the  upper  bound  on  wave  energy  as  expressed 
by  a  depth-limited  wave  height.  It  is  shown  that  in  the  shallow-water  limit 
this  leads  to  an  approximate  variation  of  H£  with  the  square  root  of  depth. 
Frequently,  the  monochromatic  limiting  value  H«j  is  used  to  provide  an  upper 
bound  on  the  wave  height  in  shallow  water.  This  report  indicates  that  such  an 
approach  can  significantly  overestimate  the  significant  wave  height.  The  tra¬ 
ditional  method  of  estimating  wave  conditions  in  shallow  water  has  been  to 
obtain  an  estimate  of  H1/3  in  some  depth  of  water,  then  refract  and  shoal 
it  into  the  shore.  At  some  point  becomes  larger  than  Hd,  in  which 

case  Hw3  Is  set  to  Hd.  This  report  Indicates,  however,  that  the  wave 
height  H£,  which  is  directly  related  to  the  wave  energy,  varies  with  h1/2 
and  is  normally  much  less  than  Hd.  Consequently,  when  the  energy  in  the  sea 
is  of  concern,  H £  should  be  used  rather  than  Hj.  If  the  maximum  individual 
wave  that  can  occur  is  of  concern  then  84  is  appropriate. 

The  method  in  this  report  also  indicates  that  the  maximum  significant  wave 
height,  H£,  in  shallow  water  in  lakes  and  bays  can  be  different  than  that  in 
the  open  ocean  because  the  cutoff  frequency,  fc,  in  the  smaller  water  bodies 
is  normally  much  higher  than  fc  for  large  ocean  storms.  Table  3  provides 
esimtates  for  as  a  function  of  h  for  an  ocean,  a  large  lake,  and  a  small 
lake  for  the  same  windspeed,  U,  of  25  meters  per  second  but  for  different 
frequencies.  Longer  waves  in  an  ocean  are  expected  to  develop  than  in  small 
lakes;  consequently,  fc  is  higher  in  the  short  fetch  cases.  The  coefficient 
a  increases  in  short  fetch  cases,  but  it  enters  H£  through  a  square  root 
relationship. 

Estimates  of  depth-limited  wave  conditions  have  traditionally  been  based 
on  linearity  of  wave  height  and  depth.  This  linear  relationship  is  well  estab¬ 
lished  for  monochromatic  waves  by  both  laboratory  and  theoretical  studies. 
Extensions  to  irregular  wave  conditions  have  relied  on  this  linear  relationship 
but  with  a  coefficient  of  about  0.4.  Figure  8  is  a  plot  of  this  variation  for 
25  October  1980  and  shows  that  in  slope  and  magnitude  this  form  is  a  poor 
predictor.  The  method  in  this  report  is  based  on  a  theory  about  spectral  shape 
and  appears  to  be  a  better  predictor.  It  should  be  noted,  however,  that 
evaluations  of  the  newer  method  must  account  for  variations  in  a  and  fc 
as  wave  conditions  change.  Hence,  simply  plotting  Kg  versus  h  or  h*/2  for 
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Squorc  Root  of  Dtpth(m,/I] 

Variation  of  significant  wavs  height,  H^,  with  the  square 
root  of  water  depth.  After  a  region  of  shoaling,  wave  height 
drops  off  linearly  with  the  square  root  of  water  depth. 
Differing  slopes  are  dus  to  variations  in  a  and  fc. 
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Figure 


7.  Estimate  of  H/  for  laboratory  conditions  (from  Seelig  and 
Broderick,  1981) .  H£  is  estimated  using  two  estimates  of 
fc,  and  a  line  H4  -  0.8  h  is  also  provided.  at  the 

toe  of  the  1:30  slope  is  13.4  centimeters  with  fp-l  -  1.47 
seconds.  A  linear  shoaling  curve  is  also  shown. 
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Figure  8.  Variation  of  wave  height  with  square  root  of  depth,  25 

October  1980,  Duck,  North  Carolina.  Solid  line  is  based 
on  measured  a  and  fc.  Dashline  represents  estimated 
band  on  monochromatic  theory  with  H  •  0.5  h. 
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Table  3.  Variation  of  H*  with  depth  for  ocean,  large 
lake,  and  small  lake  generation  cases. _ 


Depth 

Monochromatic 

Ocean 

Large  lake. 

Small  lake. 

wave,  H4  1 

*£  2 

H  3 

*1  4 

(m) 

(m) 

(■) 

(m) 

(m) 

1 

0.78 

1.5s 

1.1s 

0.9s 

3 

1  2.2 

2.6s 

1.9 

1.6 

5 

3.9 

3.4 

2.3 

1.9 

8 

6.2 

4.2 

3.0 

2.3 

10 

7.8 

4.7 

3.3 

2.6  j 

^Estimated  by  lower  limit  of  0.78  h. 

2fp  -  0.08,  fc  -  0.07,  (a/0.008’1)1/2  -  1.20. 

3fp  -  0.12,  fc  -  0.11,  (a/0. 0081) 1/2  -  1.37. 

4fp  -  0.16,  fc  -  0.14,  (o/0. 0081)1/2  -  1.44. 

sLarger  than  H4 . 

one  gage  will  show  considerable  scatter  because  of  the  tine  variation  of  a 
and  fc.  The  evaluations  of  the  method  in  this  report  have  removed  this 
constraint  by  using  a  series  of  gages  across  the  nearshore  zone. 

The  use  of  the  method  at  the  beginning  of  this  report  was  restricted  to 
spectra  of  some  breadth  such  as  storm  seas.  It  is  clear  that  nearly  mono¬ 
chromatic  waves  follow  the  linear  depth  relationship,  yet  it  is  increasingly 
clear  that  irregular  waves  do  not.  A  question  of  major  importance  not  yet 
resolved  is  how  wide  must  a  spectrum  be  before  the  waves  follow  the  relation¬ 
ships  in  this  report.  Equally  important  is  the  isolation  of  the  physics  of 
wave  motion  that  determine  these  differences.  In  a  shoaling  monochromatic 
wave,  nonlinearities  arise  which  force  the  development  of  harmonics  in  the  wave 
frequency  and  tend  to  broaden  the  spectrum,  yet  the  absence  of  other  wave 
components  may  reduce  the  transfer  energies  by  resonant  interactions.  If  the 
bottom  slope  is  sufficiently  steep,  the  evolution  of  the  swell  waves  may  be 
markedly  different  from  irregular  waves  which  may  more  easily  exchange  energy 
due  to  resonant  interactions. 


VII.  SUMMARY 

A  method  for  estimating  depth-limited  significant  wave  height,  H £, 
based  on  a  theoretical  fora  for  the  shape  of  shallow-water  storm  wave  spectra 
was  presented.  The  method  requires  an  estimate  of  the  peak  frequency  of  the 
wave  spectrum,  fp;  knowledge  of  the  Phillips'  equilibrium  coefficient,  a; 
and  water  depth,  h.  A  method  for  estimating  a  based  on  information  about 
the  peak  frequency  of  the  sea  spectrum  is  also  ^iven.  The  results  indicate 
that  the  depth-limited  significant  wave  height,  H£,  based  on  the  energy  of 
the  sea  state  is  generally  less  than  the  depth-limited  monochromatic  wave 
height,  Hq.  The  depth-limited  wave  height  defined  as  4.0(E)1/2  appears 
to  be  related  to  the  square  root  of  depth. 
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